On the variation of the gauge couplings during inflation 
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It is shown that the evolution of the (Abelian) gauge coupling during an inflationary phase of 
de Sitter type drives the growth of the two-point function of the magnetic inhomogeneities. After 
examining the constraints on the variation of the gauge coupling arising in a standard model of 
inflationary and post-inflationary evolution, magnetohydrodynamical equations are generalized to 
the case of time evolving gauge coupling. It is argued that large scale magnetic fields can be copiously 
generated. Other possible implications of the model are outlined. 
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I. INTRODUCTION 

Prior to the formation of the light elements (taking 
place at a temperature of roughly 0.1 MeV) the gauge 
coupling could have been dynamical . Two examples 
in this direction are models involving extra-dimensions 
H and string motivated scenarios 

Suppose that a minimally coupled (massive) scalar 
field (j> evolves in a conformally flat metric of Friedmann- 
Robertson- Walker type: 



ds 1 



G lJLU dx^dx v 



a 2 (rj)[drj 2 — dx 2 



where a(rj) is the scale factor, 77 the conformal time co- 
ordinate [related to cosmic time t as a{rf)drj = dt] and 
G^y the space-time metric. The field <f> is not the infla- 
ton but it evolves during different cosmological epochs 
parametrised by a different form of a(rj). Typically the 
Universe evolves from an inflationary phase of de Sitter 
(or quasi-de Sitter) type towards a radiation dominated 
phase which is finally replaced by a matter dominated 
epoch. 

The evolution equation of <fi in the background given 
by Eq. (|l|) can be written as 



<f>" + 2H4>' + m 2 a 2 (t) = 0, H 



a 
a 



(2) 



where the prime denotes derivation with respect to the 
conformal time coordinate and TL is the Hubble factor in 
conformal time related to the Hubble parameter in cos- 
mic time H — aj 'a as Ha — TL (the dot denotes derivative 
with respect to cosmic time). 

If cf> evolves during an inflationary phase of de Sitter 
type the scale factor will be 



a(rj) 



(3) 



where —771 marks the end of the inflationary phase. If 
m 2 a 2 <C TL (i.e. m <C H) during inflation, according to 
Eq. (||), <f> relaxes as <j> ~ (—77/771 ) 3 for 77 < —771. 

Suppose, as an example, that (f> is coupled to an 
(Abelian) gauge field 



(4) 



The normal modes of the hypermagnetic field fluctua- 
tions Bi(x,rj) are bi(X,rj) = 4>(j])Bi(x,rj) and their cor- 
relation function during the de Sitter phase can then be 
written as 



Gij(r,r)) = 



d 3 k 
(27js 



P ii (fc)&(fc,77)6*(fc77)e lfef , (5) 



where 



Pii = % 



hi kj 



(i) The normal modes b(k, rf) will evolve as 



b" + [k 2 -Ub = 0. 



(6) 



(7) 



Using now the fact that 4> ~ rf the correlation function, 
during the de Sitter phase grows as 

QiA^r}) ^ {-ri)-\ (8) 

for 77 — > 0~ (i.e. t — > 00). Thus, gauge field fluctuations 
grow during the Sitter stage. Furtheremore, from Eq. (^) 
the magnetic energy density pe (r, rj) [related to the trace 
of Gij(r,r])] also increases for 77 — ^ — . Consequently, 
since the magnetic energy density can be amplified dur- 
ing a de Sitter-like stage of expansion, large scale gauge 
fluctuations pushed outside of the horizon can generate 
the galactic magnetic field. 



II. EVOLUTION OF THE GAUGE COUPLING 

The only gauge coupling free to evolve, in the present 
discussion, is the one associated with the hypercharge 
field leading, after symmetry breaking, to the time vari- 
ation of the electron charge. In a relativistic plasma 
the conductivity goes, approximately, as T/a cm where 
ce cm is the fine structure constant |(|. If a om depends on 
time also the well known magnetohydrodynamical equa- 
tions (MHD) H will have to be generalized, leading, ul- 
timately, to different mechanism for the relaxation of the 
magnetic fields. 

If the evolution of the Abelian coupling is parametrised 
through the minimally coupled scalar field <f>, the pos- 
sible constraints pertaining to the evolution of 4> are 



1 



translated into constraints on the evolution of the gauge 
coupling. Massless scalars cannot exist in the Universe: 
they lead to long range forces whose effect should appear 
in sub-millimiter tests of Newton's law. Consequently, 
the scalar mass should be, at least, larger than 1CP 4 
eV otherwise it would be already excluded 0. Mas- 
sive scalars are severely constrained from cosmology ||[)| . 
When the scalar mass is comparable with the Hubble rate 
(i.e. m ~ H) the field starts oscillating coherently with 
Planckian amplitude and <f> decays too late big-bang nu- 
cleosynthesis (BBN) can be spoiled [0 . 

Initial conditions for the evolution of <f> are given during 
a de Sitter stage of expansion. Thus, the homogeneous 
evolution of d> can be written as 



77 < -T]x, 



(9) 



where <f>\ is the asymptotic value of </> which may or may 
not coincide with the minimum of V{4>); <p2 is also an 
integration constant. Without fine-tuning </>i and </>2 
both coincide with Mp. During the inflationary phase 
to <C Tii/ai = Hi where, Hi < ICP 6 Mp is the curva- 
ture scale at the end of inflation. 

After rji the Universe enters a phase of radiation domi- 
nated evolution (possibly preceded by a reheating phase) 
where the curvature scale decreases. When H m ~ m the 
scalar field starts oscillating coherently with amplitude 

01- 

During reheating the scale factor evolves as a(rj) ~ rj a 
so that, in this phase, <f> relaxes as 



7/1 



-2d 



(10) 



where r/ r marks the beginning of the radiation dominated 
phase occurring at a scale H T > m. In the case of matter- 
dominated equation of state during reheating a ~ 2. 

For r\ > rj r , the evolution of the field </> can be exactly 
solved (in cosmic time) in terms of Bessel functions 

4>{t) ~ a- 3/2 {t)Vmt[AY 1/4l (mt) + BJ 1/4 (mt)] , (11) 

where A and B are two integration constants. From Eq. 
(pd|), cf> ~ constant + rp 1 for H > m, and it oscillates 
for H < to. When H < m the coherent oscillations of 
4> start and their energy density decreases as a~ 3 . The 
curvature scale H c marks the time at which the energy 
density stored in the coherent oscillations equal the en- 
ergy density of the radiation background, namely 



H r Mp{ — ) - 777 <M — ) ' 



(12) 



where 7/ m correspond to the times at which H ~ to. From 
Eq. P) 



(13) 



where (p — (f>i/Mp and ^ = m/Mp. The phase of domi- 
nance of coherent oscillation ends with the decay of (j) at 
a scale dictated by the strength of gravitational interac- 
tions and by the mass m, namely 



Ha 



(14) 



In order not to spoil the light elements abundances we 
have to require that H^ > H ns implying that to > 10 
TeV. 

In order not wash-out the baryon asymmetry pro- 
duced at the electroweak time by overproduction of en- 
tropy |13,pJ H<f, > H cw may be imposed. Since H cw ~ 



/NeffT^/Mp [where A^ e g = 106.75 is the effective num- 
ber of (spin) degrees of freedom at T cw ~ 100 GeV] we 
obtain to > 10 5 TeV. Notice, incidentally, that the time 
variation of the gauge couplings during the electroweak 
epoch (possibly in the presence of a hypermagnetic field 
|]l3| ) has not been analyzed and it may be relevant in 
order to produce inhomogeneities at the onset of BBN 
& 

The inhomogeneous modes of <f> should also be taken 
into account since we have to check that further con- 
straints are not introduced. In order to find how many 
quanta of the field <fi are produced by passing from the 
inflationary phase to a radiation dominated phase let us 
look at the sudden approximation for the transition of 
a (77) ]l2|| . Consider the first order fluctuations of the 
field 4> 

(j)(x,r]) = <j){r]) + 5(j){x,rj), (15) 

whose evolution equation is, in Fourier space, 

ip" + 2Hip' + [k 2 + to 2 a 2 ]V> = 0, (16) 

where ip(k, rf) is the Fourier component of 84>{x, rf). 

In the limit krji <C 1 the mean number of quanta cre- 
ated by parametric amplification of vacuum fluctuations 

ira is 



n(fc)~|c_(fcr = #77in A ( Wi 



-1/2 



(17) 



where q is a numerical coefficient of the order of 10~ 2 . 
the energy density of the created (massive) quanta can 
be estimated from 



dp,p = -^^mn(k), 



(18) 



where u> = k/a is the physical momentum. In the case of 
a de Sitter phase (A = 3/2) the typical energy density of 
the produced fluctuations is 



-1/2 / x 3 
777 \ / CL\ 



Hi 



(19) 
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Also the massive fluctuations may become dominant and 
we have to make sure that they become dominant after 
(f> already decayed. Define H* as the scale at which the 
massive fluctuations become dominant with respect to 
the radiation background. The scale H* can be deter- 
mined by requiring that p^,(?y*) — p 7 (ry*) implying that 



-1/2 



Hi Mp 



which translates into 



(20) 



(21) 



where e = H±/Mp. In order to make sure that the non- 
relativistic modes will become dominant after ip already 
decayed we have to impose that < which means 
that m > 10 2 TeV for B x ~ 1(T 6 M P and which is less 
restrictive than the other constraints previously derived 
in this paper. 



III. MAGNETOGENESIS 

The full action describing the problem of the evolution 
of the gauge coupling in this simplified scenario is 



S 



d 4 x 



Using Eq. (pi) the equations of motion become 



,dV 



1 df 



B z -E 2 



dB - - 
— = -X7xE 



d_ 

drj 



J = /(0)V x B, 



V-B = 0, J = a{E + vx B) 



(22) 

(23) 
(24) 

(25) 
(26) 



(B = a 2 B, E — a 2 £; J — a j; a = <r c a; B, £, j, a c 
are the flat-space quantities whereas B, E, J, a are the 
curved-space ones; v is the bulk velocity of the plasma). 

In Eqs. (|2^) ( p6|) the effect of the conductivity has 
been included. The current density [present in Eq. ( p3] ) 
with a term (dj a /d(j>)A a ] has been eliminated by the us- 
ingMaxwell's equations. During the inflationary phase, 
for r) < —rji, the role of the conductivity shall be ne- 
glected. In this case the evolution equation for the canon- 
ical normal modes of the magnetic field can be derived 
from the curl of Eq. (p5|) with the use of Eq. (p3) : 



(27) 



where b = \fjB. For r\ > —r\\ the effect of the conductiv- 
ity is essential. Therefore, the correct equations obeyed 
by the magnetic field will be the generalization of the 
MHD equations whose derivation will be now outlined. 

MHD equations represent an effective description of 
the plasma dynamics for large length scales (compared 
to the Debye radius) and short frequencies compared to 
the plasma frequency. MHD can be derived from the ki- 
netic (Vlasov-Landau) equations and the MHD spectrum 
indeed reproduces the plasma spectrum up to the Alvfen 
frequency [§. MHD can be also derived j(| by neglecting 
the displacement currents in Eq. (Ea): 



/V x B = J + f'E. 



(28) 



By now using the Ohm law together with the Bianchi 
identity we get to 

f -> 1 
(H^-)B' = Vx(i)xB)|-V 2 B (29) 
s/ s 

which is the generalization of MHD equations to the case 
of evolving gauge coupling. The quantity s = a/ f is 
constant. The reason for this statement is the following. 
The rescaled conductivity, 



T 



(30) 



where a cm ~ f^ 1 . Therefore a/f = s with these rescal- 
ings, is constant. Taking now the Fourier transform of 
the fields appearing in Eq. (^9|) the solution, for the 
Fourier modes, will be 



B i (fc,ry)=S i (fc,? 7l )e / ' . 
Consider now, as an example, 



/(</>) 



'1x2 



(31) 



(32) 



For 7/ < — 771 the solution of the evolution equation of the 
magnetic fluctuations is, from Eq. (127 



b{k, rj) = N^/hqH^ (fcry), N = 



'kit 



-if (l+2v) 



(33) 



where TV has been chosen in such a way that b(k, rf) — > 
^/k/2 e -ikv f or v _^ _oo. Using Eq. © v = 5/2. 

For 771 < r) < r) r the Universe is reheating. During 
this phase the conductivity is not yet dominant and the 
fastest growing solution outside the horizon is given, in 
the case of Eq. ©, by 



(34) 



where we assumed, for concreteness, that a = 2 in Eq. 
@. For 77 > n r Eqs. @ should be used. 
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The typical present frequency corresponding to the end 
of the inflationary phase is given, at the present time 770, 
by 

CJxfao) ~ 1CT 4 Tdcc e^TT v ,-2/3 j (35) 

Notice that 10~ 4 T dcc = 100 GHz where Tdec is the de- 
coupling temperature. The typical frequency correspond- 
ing to the onset of the radiation dominated phase is given 
by 



Wl(lJo) 



(36) 



For 77 > rj T the conductivity dominates the evolution 
and using Eq. we can estimate the trace of the two- 
point function (n) 



pB(r,rj) = / PB(k,Tj) 



sin kr dk 
kr k ' 



with 

PB{k,Tj) 

Thus, in terms of 

r(w) 



\b(k, V )f 



Pb(u,v) 
Pi, (17) 



and using Eqs. d33|)— (J2 



fs(^G,?7o) = C(v, u>i, w r )(— ) 5 2 "^"(wg) 



(37) 



(38) 



(39) 



(40) 



where luq ~ 10 14 Hz is the present frequency corre- 
sponding to a Mpc scale and 



o2i/-2 /, , \ 4(a+l) 

7T Z V UJ T 



[fflO*?)*]^. (41) 



with w a - and T - 10~ 13 GcV. In Eq. © 

uj m and w c are, respectively, the present values of H^, 
H m and H c . Using the notation of Eq. ([35]) we have that 



2 £ -l/(a+l) ^1/2 ^(l-a)/(2a+2) Wi 
, £7/6 2/3 e -l/(a+l) ^1/2 ^(l-a)/(2a+2) ^ _ 



(42) 



All the frequencies are evaluated at the time t/q. Using 
now the previous equations, 



t b {uj g ) ~ e 4 C" 



(43) 



This result should be confronted with typical values of tb 
required in order to explain galactic (and possibly inter- 
galactic) magnetic fields. 



Large scale magnetic fields are a well known compo- 
nent of the interstellar medium. The most reliable es- 
timates available today rely on Faraday rotation of ra- 
dio signals. With these techniques the magnetic fields of 
different galaxies have been measured both within and 
beyond our local group Recently magnetic fields 

in clusters have been shown to possess a magnetic field 
which is larger than previously thought and of the or- 
der of the n Gauss [fl7| . These measurements have been 
made possible through the combined analysis of x-rays 
bright Abell clusters performed with VLA telescope and 
ROSAT satellite full sky survey. 

Through differential rotation of the primeval galaxy 
some small magnetic seeds can be substantially amplified. 
Since by the time of the formation of the galaxy the gauge 
coupling is frozen the correct equation describing large 
scale magnetic will be exactly Eq. (^9|) with /' = and 
s — > a . In this equation an instability develops when the 
dynamo term (containing the bulk velocity) dominates 
over the diffusivity term (containing the conductivity). 

Taking into account that the rotation period is of the 
order of 10 s yrs and that the age of the galaxy is of the 
order of 10 10 yrs the typical amplification which could be 
obtained through the dynamo mechanism is of the order 
of 30 e-folds, namely, 13 orders of magnitude. This im- 
plies that if today galactic magnetic fields have /1 Gauss 
strength, at the end of gravitational collapse they should 
have been as small as 10 -19 Gauss in order to turn on 
the dynamo mechanism. 

When the primeval galaxy collapses (from a typical 
scale of 1 Mpc down to a scale of 30 kpc) the frozen-in 
magnetic flux increases of roughly 4 orders of magnitude. 
This is because the mean density prior to collapse is of the 
order of the critical density, whereas, after collapse, the 
mean density of the galaxy is approximately six orders of 
magnitude larger. These estimates have been obtained 
for h = 0.65, fl\ = 0.7 and O mat t e r = 0.3. 

Taking now into account this last point the primordial 
seeds should be as strong as 10~ 23 Gauss over a typical 
scale Wq 1 prior to gravitational collapse. In terms of 
rs{(jj) this means that in order to turn on the dynamo 
mechanism we should have Il5(| 



r B (u G ) > 10 



-34 



(44) 



Various mechanisms have been proposed so far in order 
to explain the magnetic field of the galaxy |18|. Recently 
|]l9| , it was pointed out that the evolution of gauge cou- 
plings in a Kaluza-Klein context can offer interesting pos- 
sibilities for the generation of primordial magnetic fields. 
The present analysis could be viewed as a realization of 
that proposal in a geometry compatible with a de Sitter 
stage of inflation. 

In the case of clusters the dynamo mechanism is more 
problematic. Indeed on one hand clusters rotate less than 
galaxies and, on the other hand, the mean density of mat- 
ter is smaller. Therefore, larger values of rs are required 
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FIG. 1. Magnetogenesis requirements are illustrated. The 
thin line b corresponds to rB(wG) ~ 10 -34 . The dashed line 
a corresponds re(tJG) ~ 10 -20 . The (horizontal) thick line 
correspond to e = 10 -6 . The (diagonal) thick line corresponds 
to e ~ (. Magnetogenesis is possible when the parameters lie 
in the shaded area. 



in order to successfully implement the dynamo mecha- 
nism for clusters. 

Using Eqs. (^) and (|44|) the allowed region in the 
parameter space of the model can be obtained by tak- 
ing into account the constraints discussed in the previous 
section. In Fig. |l| the shaded area illustrates the region 
where magnetogenesis is possible. To be consistent with 
inflationary production of scalar and tensor fluctuations 
of the geometry e < 10~ 6 should be imposed. Thus, in 
Fig. [l] the parameters should all lie below the (horizon- 
tal) thick line. Moreover, since H T < Hi, £ < e. Recall 
that C > 10 -15 in order not to affect the nucleosynthesis 
epoch [see the dot-dashed line in Fig. Q. This require- 
ment comes about since H x > m and £ > 10 . 

In this paper the possible phenomenological implica- 
tions of the evolution of the gauge coupling have been 
analyzed in the case of a specific model of background 
evolution. It has been shown that if the gauge coupling 
is related to a massive scalar field minimally coupled to 
the geometry the phenomenological constraints related 
both to the evolution of the massive scalar and to the 
evolution of the gauge coupling can be satisfied. In a 
specific example the large scale magnetic fields produced 
with this mechanism have been computed. It has been 
shown that they can be large enough to seed the galac- 
tic dynamo mechanism. They can be also relevant for 
the origin of magnetic fields in clusters. Further work on 
these possibilities is in progress. 
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